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r*j . Abstract 

We consider rectangular Wilson loops in certain non-equilibrium quantum 
states in M = 4 SYM at weak coupling, prepared with a quantum quench. 
We find that in the ladder approximation, the Bethe-Salpeter equation can be 
reduced to solving a massive 1+1 dimensional wave-equation with a leaking 
boundary condition leading to a quasinormal behavior analogous to what is 
found in studying dynamics of fields in black hole backrounds. Furthermore, 
we find that the Wilson loops with size L thermalize after time T ~ L/2 in 
qualitative agreement with gravitational calculations in falling shell models of 
thermalization. 
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1 Introduction 

There is a large amount of evidence that M = 4 super Yang Mills theory is secretly a 
string theory in AdS xS 5 ,a quantum theory of gravity. One of the main differences 
between gravitational theories and non-gravitational theories is the amazing property 
that gravity can dynamically alter the causal structure of the theory. A manifestation 
of this difference is black holes, which have the property of absorbing field fluctuations 
never letting them back to the outside spacetimqj. 

Static black hole geometries are dual to thermal states in the field theory. A 
qualitative understanding of many of the features of black holes can been understood 
in terms of large N gauge theory in a thermal state. The purpose of this work is to 
take a small step towards understanding black hole formation from the perspective of 
M = 4 SYM (for other work in related directions see [IHS])- We consider deforming 

1 This is when considering large black holes in AdS which do not evaporate. 
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SYM with a time dependent source for the Konishi (mass) operator Tr^ 1 ^ 1 . More 
specifically, we start with a non-vanishing value of the source and suddenly turn it off. 
This way the subsequent time evolution is generated by the usual SYM Hamiltonian. 
The state is a high energy state which can be expected to thermalize. 

In the dual string theory picture, fast varying time dependent sources are expected 
lead to the creation of a shell of energy density near the boundary, which subsequently 
falls into AdS forming a black hole [HIE]- In our case the operator corresponds to 
a massive string mode, which is beyond the above considerations. Still one might 
expect a similar qualitative picture of black hole formation to apply. 

Correlation functions [8j[9], entanglement entropy [TO] and Wilson loops [8j|9] have 
been studied in such falling shell backgrounds. The behavior of correlation functions 
and entanglement entropy in conformal (and other) field theories subject to a quantum 
quench (a time dependent coupling) have been studied for example in [TTHT3] . In this 
work we study Wilson loops in the non-equilibrium state described above. 

Perhaps the main qualitative lessons arising from the collapsing shell calculations 
is that correlations separated by a spacelike distance L thermalize maximally fast in 
time T f» L/2 [HIE], with a memory on the initial state being of the order of the time 
scale of the lowest quasinormal mode [TTt[T5] . This is the same time scale as one finds 
in considering linear response around thermal equilibrium. 

To calculate the Wilson loop we use two approximations. The first one was in- 
troduced in [16J to isolate the scalar ladder Feynman diagram contributions to the 
Wilson loop. We will call this the ladder approximation. Somewhat surprisingly the 
results from the ladder approximation can be continued to strong coupling (when 
the coupling constant is properly defined through analytic continuation) and so far 
impressive agreement with string theory results is found [T6T4T9] . Also first order 
corrections to the leading ladder approximation are found to agree with the strong 
coupling results [20] . 

As a second approximation we take a limit in which the initial value of the scalar 
mass (the source of the Konishi operator) is large as compared to the inverse sep- 
aration 1/L of the timelike segments in the Wilson loop. This approximation was 
introduced in [32JQ3] to study scalar two point correlations in such a non equilibrium 
state. Following the terminology of [12l[T3] we wm C& U this the deep quench approx- 
imation. In this approximation the two point function of the scalar field simplifies 
drastically allowing for analytic calculations. 

By summing up the ladder diagrams we end up with a Bethe-Salpeter equation 
that can be written as a massive 1+1 dimensional wave-equation with a "leaking" 
boundary condition, at an emergent boundary, that lets field fluctuations cross the 
boundary and never come back. This leads to a quasinormal behavior of the Wilson 
loop. This is qualitatively quite similar to the quasinormal behavior of fields in a 
black hole background. We are not aware of calculations of rectangular real time 
Wilson loops in black hole backgrounds, but heuristically it seems plausible that the 
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Wilson loops decay exponentially as a function of the temporal length. 

The leaking boundary condition inevitably leads to a fast thermalization of the 
Wilson loop in the same way as fluctuations of fields in a black hole background decay 
fast. This can be contrasted to expectations in large N vector models that are found 
to thermalize through a somewhat slow time evolution with several well separated 
time scales (see e.g. section 4.1.4 of and references therein). This is in qualitative 
agreement with the recent results on duals of large N vector models which suggest 
that the thermal ensemble is not described by a black hole state [22], but rather by 
a thermal gas, so that in the bulk thermalization is expected to happen far slower 
through interactions in the bulk gas of excitations. 



2 Setting up the non-equilibrium state 

We start from M = 4 SYM with a time dependent source turned on for the Konishi 
operator Tr($ / $ / ), 



S = -^ J d 4 xTr(^F 2 + (D&) 2 +m 2 (t)(<S> 1 ) 2 - $ J ] 2 + fermions) . 



(1) 



A convenient choice for calculations is to choose m 2 (t) = 9{—t)m^. We believe that a 
more generic choice of m{t) interpolating between mo and will lead to qualitatively 
similar results. 

We will begin from the ground state at times t < 0. In the g — > limit this state 
is annihilated by the corresponding annihilation operators of the free SYM fields. In 
momentum space, and in the Heisenberg picture, the free equations of motion for $ 
at times t < are solved by 

k) = 9 (Be- iuJo{k)t + B ] e iu)o{k A . (2) 

After t > 0, the corresponding solution is 

k) = - 9 ( Ae~ iuj{k)t + A^e luj(k)t ) . (3) 

The factors of the coupling constant g are chosen for normalization. Above we have 
introduced the notation uo(k) = \Jk 2 + ml and u(k) = \k\. Matching these two gives 
us a Bogoliubov transformation relating the two sets of creation and annihilation 
operators 

A = aB + (5B\ A ] = a*B ] + (5*B (4) 
where 
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Thus, the quench induces a non-vanishing particle number 



(N k ) = |/3| 2 = ^y^^f • (6) 
Ak^k 2 + ml 

This can be compared to a thermal ensemble (Nk) = [eP k — l) -1 - At small momenta 
these have the same form (Nk) ~ l/(kj3), if we would identify m /4 as an effective 
temperature. The particle number differs significantly from the thermal one at large 
momenta as (Nk) ~ mo/(16A; 4 ), which in the thermal case would vanish as e~^ h . 

The two point function in this time dependent state can be worked out by using 
the fact that the quench state satisfies B\ip) = 0, which gives 

{& a (t,x)$ J b (t?,0)) = g 2 S ab S IJ x, 



d 3 k e lk ' x 
(2tt) 3 2\k\ 



a \2 e -i\k\{t-«) + \p\2 e i\k\{t-t>) + a(3e - im+ t>) + a *p* e i\k\(t+t>)\ s (?) 



where we have written down the SO (6) and the color indices explicitly. In what 
follows we will again suppress these indices. In the rest of this paper we consider a 
large quench limit (as was used in [13J), where mo 3> 1/L, where L is a characteristic 
length scale that we are interested in. More precisely, when \x\ — (t — t') ^> l/m , 
the propagator can be simplified 

($(t,x)$(t / ,0)) w g 2 I™ -^—sm(k\x\)^( cos kit - t') - cos kit + if)). 

J 87i z \x\ k V / 

By performing the integral using the identities in Appendix |X] we obtain 

K(t,x;t',0) = 0)> « ^i-e(N, t,t'), (8) 

1D7T \X\ 

where we have defined 
0(\x\,t,t') = -(sgn(\x\+t-t')+sgn(\x\-t+t')-sgji(\x\+t+t')-8gn(\x\-t-t')). (9) 

This function is visualized in Figure [TJ The two point function (JS} is the Wightman 
function, but it should be noted that in the deep quench approximation, the expec- 
tation value of the field commutator is subleading in l/m , making the Wightman 
functions and time ordered two point function equivalent, in the leading order. 

The approximation leading to OH]) breaks down near the light cone \x\ — (t — 
t') = 0. Instead, the two point function has the same singularities on the light cone 
as in the vacuum state. These singularities are likely irrelevant for the qualitative 
behavior of the Wilson loop. We do not have methods to deal with the singular parts 
quantitatively. Thus, we will ignore them for most of the calculations and comment 
on their effects in Appendix O 



5 




Figure 1: The function Q(\x\,t,t r ) visualized on t,if plane for a fixed value of \x\. 

3 Climbing the first ladder 
3.1 Specifying the Wilson loop 

To specify a real time Wilson loop one needs to specify a curve in spacetime. In the 
current paper we work in Minkowski spacetime R 3 ' 1 . For most of this paper we will 
choose to study rectangular Wilson loops whose sides have a temporal length T and a 
spatial length L, shown in Figure [21 In this paper we consider the Maldacena- Wilson 



where P denotes path ordering of the matrices. This is the Wilson loop that sources 
a fundamental string in the dual string theory [23]. To define the Wilson loop ffTUj) 
one has to specify a direction of the symmetry breaking VEV of the (external) W- 
boson in the R-symmetry group SO (6). We will choose the relative angle between the 
upper timelike segment and the lower timelike segment to be <p. While the relative 
angle between the spacelike segments and the lower segment is chosen to be (p/2. The 
choices of the SO (6) orientations are displayed in Figure |2j 

In what follows we use the approximation introduced in [16] which isolates the 
scalar ladder diagrams as the leading contribution to the Wilson loop. For n $ 
propagators going between the timelike segments we have a factor of (cos(ip)) n . For 
the planar ladder diagrams there is also a corresponding factor of A n . More general 
diagrams with more interactions on the internal lines have more powers of A. This 
suggests a limit which picks up the scalar ladder diagrams 



loop [23] 




N 



(10) 



<p — ii?, d — > oo 



(11) 
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while keeping 

A = \e»\ (12) 

fixed. Diagrams with lines running between spacelike and timelike lines are suppressed 
in this limit as Ae -1?//2 , while the lines running between the spacelike segments are 
suppressed as Xe~^. All the corrections to the free scalar propagator are suppressed 
by powers of Ae _1? . It is likely that the approximation of neglecting such diagrams 
will not commute with the large distance limit. 



(l,o) e 2 (l,t) e 2 




Figure 2: The rectangular Wilson loop, with SO(6) orientations denoted as 0's. 



There is some freedom in defining a real time Wilson loop. One can consider 
several different orderings of the operators (see e.g. (SUES]). One could for example 
consider a time ordered operator ordering, or a path ordered one with some choice 
for the starting point of the path. In the leading order deep quench approximation, 
the time ordered two point function is equivalent to the Wightman functions, making 
these choices to give the same result for the Wilson loop. Thus, we can simply ignore 
this subtlety in this work. 

3.2 1-loop warmup 

First consider the first order correction to the Wilson loop as a warmup. This is given 
by 

re> « i - e ^ Tl ^ aTb) J\ tl £ dh^^mx^ix,)}. (is) 

Using T a T a N/2 and the two point function (jSJ) we get 

f dt [ ds^e(L,t,s). (14) 
16vr Jo Jo L 

Rather than giving the full result for the integral we note the most important things 
about the answer. Firstly the whole correction vanishes for IT —L < 0, since Q(L, t, s) 
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vanishes within the corresponding integration region. Secondly, for large times T 3> 
L/2 we find 

(W(C)} « 1 - ^T. (15) 

Since the two point Wightman function (jSJ) is real, we immediately see that the Wilson 
loop will be real. We will come back to this point later on. 



4 Summing up the ladders 
4.1 The Bethe-Salpeter equation 

It is well known that ladder diagrams can be resummed using a Bethe-Salpeter equa- 
tion [H1HB1I25]. 



t 




s 



Figure 3: The Bethe-Salpeter equation resumming the ladder diagrams. 



The Bethe-Salpeter equation in question is graphically shown in Figure [31 As an 
equation it reads 

n r s r T 

r(S,T) = l-—cos(p ds dtK(t,L;s,0)F(s,t), (16) 
2 Jo Jo 

where K is the two point function (jBD. Taking two derivatives of (TTSj) gives a differ- 
ential equation for V 

d { sdT ' = ~ cos (pK(T, L- S, 0)T(S, T). (17) 
This equation should be solved with the boundary conditions 

r(s,o) = i, r(o,T) = i, (is) 

These boundary conditions follow from the integral form ( TloT) . 

A more familiar form of the equation is obtained by defining x = S — T and 
y = S + T. This way the equation becomes 

(d 2 y -d 2 x +m 2 eff )r(x,y) = 0, (19) 



L 



-L 




Figure 4: The Klein-Gordon equation (1191) is massless in the regions I and III. The 
effective mass is non-vanishing in the region II. The dashed lines correspond to the 
light-cone where the initial conditions are specified. 



which is nothing but the Klein-Gordon equation with a spacetime dependent mass 
term 

The form of the mass term is shown in Figure HI The boundary condition for the 
Klein-Gordon equation is given at the light-cone x = ±y with y > 0. 

4.2 The horizon effect 

First consider solving the BS equation ( Tl9l) in the region I shown in Figure HI In this 
region, it reduces to a massless Klein-Gordon equation. 

The massless Klein-Gordon equation is most conveniently analysed using the orig- 
inal coordinates S and T. The equation reads 

d s d T T = 0. (21) 

The general solution is simply 

T = f(T) + g(S), (22) 

where the functions / and g are fixed by the boundary conditions T{S, 0) = T(T, 0) = 
1, which gives 

f(P) + g(T) = l, f(S)+g(0) = l. (23) 
Thus, / and g must be constants, that satisfy / + g = 1 , which fixes 

T(S,T) = f + g = 1. (24) 
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This immediately tells us that in the early time region y < L, where m e ff = 0, the 
solution is simply T = 1. 

This means in particular that all of the ladder diagram corrections to the Wilson 
loop vanish when T < L/2. This is analogous to the "horizon effect" noticed for 
correlation functions in [12] and for entanglement entropy in [llj. The "horizon 
effect" reflects the fact the in order for the quench to generate correlations between 
distant points, one at least has to wait for the time T = L/2 for the correlations 
between spacelike separated operators to be causally affected by the quench. 



4.3 A Heuristic picture of the decay of the Wilson loop 

To understand the qualitative behavior of the solutions of the Klein-Gordon equation 
(11 9p for later times, it is convenient to introduce a 1 + 1 dimensional energy momentum 
tensor for the field T 

t„ v = -^u(dr) 2 + 9 M ra„r - ^ % „m 2 eff r 2 . (25) 

The four momentum density T M o is conserved separately for times y < L and y > L, 
while its value jumps at y = L. For y < L the solution that satisfies the appropriate 
boundary conditions is T = 1. Thus, the energy of the fields vanishes for y < L. 
As the time passes y = L, the mass term gets a non-vanishing value and the system 
obtains an energy 

E = J T m = ^ 2 L. (26) 

This will make the field oscillate on top of the potential hill in the region — L < x < L 
where the mass is non-vanishing. Ignoring the regions x > L and x < —L, this 
oscillation would not be decaying but would have the form T = cos(fi(y — L)) . Taking 
into account the fact that the field can leak out of the region — L < x < L leads to 
a decay of the solution. To get a more quantitative feel of how the field behaves we 
make an ansatz that the field inside the massive region has the form 

A(y)cos(fi(y-L)), (27) 

where A(y) is assumed to be a slower varying function than the cosine. Taking into 
account the fact that the energy in the field oscillations can escape from the massive 
region we obtain an outgoing energy flux 

T yx = d x Td y T oc fi 2 A 2 , (28) 

where we used the fact that the waves outside — L < x < L satisfy a massless wave 
equation so that d x T = ±d y T at x = ±L. This decreases the total energy inside the 
massive region El oc fi 2 A 2 L with a rate 

^ oc -^ 2 A 2 , (29) 
dy 



10 



which gives a differential equation for A(y) that is solved by 

A(y)oce- c ^. (30) 
This leads to the following form for the Wilson loop 

(W(C)) oc cos(2 y uT)e- cT/i . (31) 

The coefficient c is a constant factor (meaning independent of T) that is not deter- 
mined by our heuristic calculation. Thus, we see that the qualitative behavior of the 
Klein-Gordon field is the following. There are oscillations in the field in the region 
—L < x < L that escape through the boundaries x = ±L leading in to a loss of energy 
and an exponential decay of T. This behavior is characteristic to a system exhibiting 
quasinormal modes. In the next section we make this connection more precise. 

At first one might be surprised that the real time Wilson loop does not have the 
usual e~ tTE behavior, but rather behaves as e~ Tl ' . Such an exponential decay is often 
a sign of an instability [24"||25] . In the current case it could be interpreted as signalling 
that mesons become unstable in an energetic bath of particles as was found in thermal 
equilibrium in pure Yang-Mills theory in [25J. 



4.4 Reducing the Wilson loop to a quasinormal mode prob- 
lem 

To obtain a fully quantitative solution of f JT9|) it is convenient to solve the fields outside 
the potential hill —L < x < L. This can be done by noting that the the solution in 
the region x > L is 

r = T + = g + (y-x) + l-g+(0), (32) 
and correspondingly in x < — L 

r = r_ = /_(3/ + x) + i-/_(o). (33) 

The only data in these regions is the amount of outgoing radiation going towards x = 
±oo. Recalling that the Wilson loop is obtained as the value of T at x = 0, suggests 
we do not need to know the form of the outgoing radiation and can concentrate on 
the region — L < x < L. This can be done by noting that at x — ±L, F satisfies 

(d x ± d y )T = 0. (34) 

When y > L, the effective mass term is independent of y, which allows us to perform 
a Fourier transform 

T(y,x) = [ du)e~ iv "f(u,x). (35) 
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This way the field equation (fT9j) and the condition flMl) can be reduced into a massive 
wave equation within the region x G (0, L) together with a mixed boundary condition 

{-dl - tu 2 + fi 2 )f = 0, ^f(x = 0) = 0, a x f(x = L) = iuf (x = L), (36) 

where we used the symmetry of the problem under x — > —x. Of course equations ( |36l) 
are accompanied with the initial data that T{x,y = L) = 1. We can look for mode 
solutions to the equation (136]) which have the form 

T = A cos( v / w 2 - H 2 x) . (37) 

Imposing the boundary condition at the leaking wall x = L leads to a condition for 
the frequency u given by 

tan(v^^L) = -i JL_ (38) 

V u _ A* 

There are no real eigenvalues satisfying (|38l) and we have been unable to find solutions 
on the upper half complex u plane, as expected on physical grounds. Thus, we indeed 
find that we are dealing with a quasinormal mode problem. Generically solutions to 
fl38|) can be found numerically. The real and imaginary parts of the quasinormal 
frequencies are shown in Figure [5) The quasinormal modes exhibit a "phase transi- 

Re(w.) 



Im(w,) 1.4 
0.0 L , „ 




Figure 5: Solutions to the eigenvalue equation obtained numerically as functions of 
L. Left: The imaginary part of the lowest quasinormal frequency. Right: The real 
part of the lowest quasinormal frequency. The units in the figure are chosen so that 
Am /167r = 1. 



tion" at L = L c w 22/Amo- When L < L c , the modes are purely imaginary. While 
for L > L c , the quasinormal frequency obtains a real part leading to an oscillation 
of the Wilson loop. This follows from the eigenvalue equation (I3"8"|) . as two of the 
lowest quasinormal frequencies on the imaginary axis get closer to each other as L is 
increased, and at L = L c they collide and move off the imaginary axis. 
The lowest quasinormal mode at large L can be found analytically as 

47T 3 1 

Re(a;*) w //, Im(a;*) « — — . (39) 

Am Li 
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So at large L we find that the Wilson loop behaves as 

(W(C)) « cos(2/iT)e" T ^^. (40) 

For small L one can solve the Klein-Gordon equation ( 1T91) exactly as described in 
Appendix [Bj which leads to 

(W(C)} = e-^ T -^ T ~i\ (41) 

in agreement with the numerical solution of the eigenvalue equation shown in Figure 
[5]and also with the 1-loop result ( H~5l) . 

4.5 Numerical solutions 

One can also find full numerical solutions of the wave equation ffT9j) . This way we 
can check that the approximate solutions discussed in the previous sections indeed 
corresponds to the full real time evolution. 

The numerical solutions are obtained using Mathematica's NDSolve. For Lfi of 
order 1 we indeed find good agreement with the numerical solutions and the previous 
approximate solutions^ In particular we find that the exponential decay of T starts 
quickly after the time y = L. Heuristically this is indeed expected as the walls of the 
system leak right from the beginning, leading to inevitable exponential decay. 



KW>| 




Figure 6: An example numerical solution of (jl~9]) shown in loglog scale. The dashed 
line corresponds to the imaginary part of the quasinormal mode obtained as a solution 
to (|38l) . The figure corresponds to the choice of parameters Litlq = 167r/A 

An example numerical solution together with the exponential decay predicted by 
solving the eigenvalue equation ( )38|) are shown in Figure [6j The solution is seen to 

2 For L[i ;§> 1 our numerical methods of solving (fT9|) become unreliable. 
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decay exponentially with a rate given by the imaginary part of the lowest quasinormal 
mode. Also the oscillation period is found to be given by the real part of the lowest 
quasinormal mode to a good accuracy. 



4.6 Comparison to thermal equilibrium 

In this section we compare the non-equilibrium Wilson loops to real time thermal 
Wilson loopsjfl We will continue to work in the large $ limit and also take a high- 
temperature limit. The finite temperature real time Wightman function with tem- 
perature is given by 

K(t,x;t',0)^-f—Qp(x,t), (42) 

4p7T|X| 

where we have ignored terms that grow slower than 1//3, as (3 — > 0, as we are interested 
in comparing to the large quench results which hold in a limit of a large energy density. 
Above we have defined 

Qp(x, t) = sgn(x + t-t') + sgn(x -t + t'). (43) 



As the large temperature two point function (1421) has a very similar form as the 
non-equilibrium two point function, the subsequent steps leading to a Bethe-Salpeter 
equation are identical as before. We are lead to a Klein-Gordon like equation for the 
finite temperature Wilson loop Tp 

{d 2 y -d 2 x + m^Tf, = 0, (44) 

where now 

ml = ^Qfax). (45) 

Furthermore, all our results on exponential decay of the Wilson loops apply to the 
thermal case, with the only difference being the initial conditions, which basically 
just means that the "horizon effect" is absent from the thermal state as the thermal 
state is time translationally invariant. 

Thus, the non-equilibrium Wilson loops decay with the same quasinormal spec- 
trum as in the thermal case if we identify the temperature 1/(3 = mo/4, which means 
that they are practically thermal once the time has passed T = L/2. 



5 Spacelike Wilson loops 

In this section we consider spacelike rectangular Wilson loops, with sides x± G (0, L) 
and x 2 G (0,Z). We choose the 5*0(6) orientations as described in section 137X1 with 

3 Euclidcan time thermal Wilson loops have been studied in the ladder approximation in [27], 
with and without thermal screening. 



14 



the replacement T — > Z, so that only the ladder diagrams between the constant X\ 
lines contribute. In this case there are no ambiguities with real time time ordering as 
the loop is located at constant time. The two point functions that contribute to the 
ladders are equal time two point functions. The Bethe-Salpeter equation in this case 
becomes 



T(Z U Z 2 ) = l + ±-°9(t - L/2) / d Zl / dz 
16vr J Jo 



y/L* + { Zl -Z 2 )2' 



(46) 



To find out the time dependence of the spacelike Wilson loop, we do not even have to 
solve (l46p . We can simply note that for times t < L/2 the integral term vanishes and 
we obtain T = 1. This again gives us the horizon effect for spacelike Wilson loops. 
For later times t > L/2, the Bethe-Salpeter equation ([46]) is independent of time and 
thus T stays constant. Thus, as far as the time dependence is concerned the Wilson 
loop has the form 

(W) = e(t-L/2)f(L). (47) 

The L dependence can be worked out for example by using the methods of [17] . The 
function f(L) is again the same as one would obtain at finite temperature with the 
identification = mo/4. 
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6 Discussion 




Figure 7: The qualitative form of the results for the timelike Wilson loop. Here 
= uj r + iuji are the quasinormal frequencies, which are functions of L. 



The results of our calculations are the following. The Wilson loops exhibit a " hori- 
zon effect" similar to that found for two point functions and entanglement entropy in 
1+1 conformal field theories [TT|[T3]. More precisely this means that the Wilson loop 
is simply 1 in the region T < L/2. The "horizon effect" reflects the fact the in order 
for the quench to generate correlations between distant points, one at least has to 
wait for the time T = L/2 for the correlations between spacelike separated operators 
to be causally affected by the quench. 

Furthermore we find that the timelike Wilson loop decays with quasinormal modes 
identical to those in thermal equilibrium, with the temperature identified as T = 
mo/4. The spacelike Wilson loop is found to thermalize immediately after t > L/2. 
It should be noted that the comparison to a thermal Wilson loop presented in section 
14.61 is quite naive in the sense that we are using the free thermal two point function 
for the scalar field (and of course free two point functions in the quench state). The 
form of the two point function is known to change qualitatively at large times and 
distances. In our case the one loop corrections would generate a distance/time scale of 
the order e^/moVX which can assumed to be large in the "ladder approximation" 
$ — > oo, and also at weak coupling. For the timelike Wilson loop this is not necessary 
a big problem as at such times the Wilson loop has already practically decayed to zero. 
The spacelike Wilson loop is more sensitive to the time dependence of the two point 
function, which means that its thermalization closely follows the thermalization of the 
two point function. Thus, our calculation cannot say much about the thermalization 
of the spacelike Wilson loop when finite $ effects are taken into account. 
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The thermalization picture we have obtained is in qualitative agreement with 
string theory results of |SJ|2], where it is found that short spacelike Wilson loops 
thermalize faster than long Wilson loops, with the thermalization time approximately 
satisfying T « L/2. This is easily understood, as a long string can pass through the 
falling shell surface seeing the time dependent geometry for longer times. 

The exponential decay of the Wilson loop might seem surprising at first as com- 
pared to the usual functional form e~ lET . Such exponential decay has been noticed 
for Wilson loops for example in [22J[25] • It seems likely that in our case it is signalling 
an instability of heavy mesons in the energetic bath of particles produced by the 
quench. 

It seems plausible that real time Wilson loops indeed decay exponentially in time in 
a black hole background. It is well know that particle propagators decay exponentially 
(with quasinormal behavior) in black hole backgrounds as propability falls through 
the horizon. The Wilson loop may be thought of as an amplitude for a flat string 
lying at z = e and t = extended in x direction with a length L to be found at the 
same position at a later time t = T, keeping the endpoints fixed at the boundary. 
As long as the string is sufficiently long one would expect to find a similar leaking of 
propability through the horizon, leading to exponential decay. 
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A Some useful integrals 

To perform Fourier transforms 




(48) 



we use the following representation of theta function 




(49) 



This allows us to write 




(50) 
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Thus we have 



i n+1 



«/.<«) = — t (51) 



We will need /„ for n — 1 in which cases we get 

/i(a) = — log (a + ie) + const. (52) 

A further useful identity is 

fi(a) — fi(—a) = log(— a + ie) — log (a + ie) = msgn(a). (53) 

B Short distance approximation 

In the limit L — > we can approximate the effective mas term with a delta function 



m 2 eff = £-Q(L,x,y)^Kx0(y-L)6(x), (54) 

where we have introduced the mass scale a = \m /8ir. Denoting the region x > 
by + and x < by -, we can find a general solution to the massless Klein-Gordon 
equation 

r = f±(y + x)+g ± (y-x). (55) 

Imposing the boundary condition F(y,y) = 1 for y > L and T(y, —y) = 1 for y > L 
leads to the solution 

T = T + =g + (y-x) + l-g(0), (56) 

on the region + and 

r = r_ = /_(i/ + x) + i-/_(o), (57) 

on the region — . Requiring continuity at x = leads to 

9+{y) = f-(y) + const. (58) 

The equation of motion implies that the first x derivative of T has a discontinuity at 
x = given by 

^r(y, a: = -e) - d x r(y, a; = e) = -oT(x = 0,y). (59) 

Using d x F(y,x = —e) = d y f-(y) and (9 x r(y,x = e) = —d y f-(y) gives a differential 
equation for f-(y) 

d y f-(y) = ~f-(y), (60) 

which is solved by 

/_(</) = ioe-^ /2 . (61) 
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Imposing the initial conditions leads to 

r = e -f%-L) (s/ -L)_ (62) 

Thus, the Wilson loop in the large L limit is given by 

(W(C)) = T(y = T, x = 0) = e-^Cr-TXT-*). (63) 

Expanding the exponential in powers of T we find that the first order indeed repro- 
duces the perturbative result (TT51) . 



C Comments on the light-cone singularities 

When calculating the Wightman function ((§]) we have taken the large mo limit in a 
rather naive way, keeping only terms of the order mo in the two point function while 
neglecting the order m® contributions. At large distances/times these contributions 
can be neglected as they provide only subleading contributions to the small k singu- 
larities of the momentum space two point function. There is still a situation where 
the order m[] contributions (or even the order TOq 1 ) cannot be neglected. This is 
near the light cone singularities. There the two point function takes the same short 
distance form as in the ground state. The contributions of the near lightcone region 
to the two point function have the form 

5K^frc6 mo (\x\-t), (64) 
\x\ 

where 5 mo is a linear combination of functions peaked around a region of the size 
l/m and delta functions. These singularities induce a correction term 5K to the 
Bethe-Salpeter kernel. The Klein-Gordon equation changes by the addition of walls 
of thickness I/tuq to the edges of the potential hill at x — ±L. Still, the waves can 
easily leak out of the potential hill leading to exponential decay of the Wilson loop as 
discussed in section |4"73"1 To see more explicitly that the waves indeed pass through the 
walls we can consider a WKB approximation to the Klein-Gordon equation. Consider 
an ansatz T = e -^- iS ( x )_ T n the WKB approximation, the tunneling through the 
walls is supressed through the imaginary part 

|Im(5)| = J ° dx^m 2 eff (x) -lu 2 , (65) 

where L — xq is the classical turning point with x of the order l/m . The short 
distance form of vr? e ^ near the light cone is rr? e ^ 00 ^/{ x ~~ L), which tells us that 
the WKB integral fl65|) is finite even though it passes through a singularity, and the 
tunneling is possible and not suppressed by any large negative exponent. 
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If the Klein-Gordon equation has bound states at the singular parts of 5K, they 
can keep the value of the Wilson loop bounded from below by a factor of the order 
e~ cm ° L , for some constant c. Studying these effects is outside the scope of this paper 
as they are highly suppressed in the deep quench limit. 
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